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Abstract 



Isothermal magnetization curves are plotted as the magnetization times the 
magnetic induction, AttM ■ B, versus the applied field, H. We show here 
that this new curve is the average kinetic energy of the superconducting state 
versus the applied field, for type-II superconductors with a high Ginzburg- 
Landau parameter k. The maximum of AttM ■ B occurs at a field, H*, di- 
rectly related to the upper critical field, H C 2, suggesting that H C 2(T) may be 
extracted from such plots even in cases when it is too high for direct measure- 
ment. We obtain these plots both theoretically, from the Ginzburg-Landau 
theory, and experimentally, using a Niobium sample with T c = 8.5 K, and 
compare them. 
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I. INTRODUCTION 



The average kinetic energy in Classical Physics is a useful quantity that helps to obtain 
important information about several many-body systems, such as the Sun, whose internal 
temperature can be estimated in this way [1J . In cases when the interaction among particles 
is only partly known, and the Virial Theorem applies, the total energy can be directly 
related to the average kinetic energy, and the average kinetic energy to the equipartion 
theorem, leading to the prediction of thermal properties of the classical system. In this paper 
we demonstrate, both theoretically and experimentally, that the average kinetic energy of 
the superconducting state can be directly read from isothermal magnetization curves, for 
superconductors such that the Ginzburg-Landau parameter k is larger than a few units. 
Thus this average only involves the paired electrons that form the condensate. 

It is often the case in Physics that just expressing data in a different way is sufficient 
to bring a new insight into the problem under investigation. This is exactly the present 
situation and we claim here that an interesting and new insight is obtained if isothermal 
magnetization (4nM) data is plotted differently, namely, as the product 4tcM ■ B, where 
the magnetic induction is B = H + 4ttM and H is the equilibrium applied field. In this 
paper we show that this new curve 4ttM ■ B directly determines the average kinetic energy 
for superconductors with k > k c , with an accuracy of less than one percent for k c « 3. 

A remarkable property of the product 4ttM ■ B is that it may have several local minima 
along the mixed state. This property follows from the fact that 4ttM ■ B vanishes at both 
critical fields, H c \ and H C 2. At H&, M reaches a minimum and B vanishes, whereas at H C 2, 
B is maximum and M vanishes. Since 4irM ■ B vanishes at two fields, the two extremes 
of the mixed state, it must necessarily have an absolute minimum in between, and possibly 
several local minima. The connection between the average kinetic energy and the isothermal 
magnetization follows from the Virial Theorem for superconductivity || which states that, 

^ = Fkin + 2 F field, (1) 

where F^n and F 'field = (h 2 /87r), are the average kinetic energy and the average magnetic 
energy, respectively. The local magnetic field inside the superconductor is h and (• • •) 
represents an average value taken over the volume of the superconductor. This theorem has 
been discovered rather late for the Ginzburg-Landau (GL) theory and it has been proven 
very useful since then, because it avoids taking the derivative of the free energy for the 
determination of H. The theorem has been extended to the microscopic theory in the context 
of the quasi-classical Eilenberger theory ||. There it was found that the above theorem holds 
for temperatures well below T c and also in presence of non-magnetic impurities. Thus one 
expects that even in presence of some types of pinning the above relationship still holds. 

In this paper we obtain 4irM-B curves both theoretically and experimentally for isotropic 
superconductors. We only discuss here examples such that the symmetry along the direction 
defined by the applied field turns the theoretical problem into a two-dimensional one. In this 
case the vector notation is no longer necessary and so, is dropped hereafter. The 4irM ■ B 
theoretical curve is derived in three situations of the macroscopic GL theory. We obtain 
analytical expressions for it in the Abrikosov and the London limits, we also calculate this 
curve numerically, using the iterative scheme proposed by E.H. Brandt |3J], a very accurate 
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method that gives the periodic solution of the GL theory for any field in the mixed state. 
Using the iterative method we show that for k larger than a few units the curve AtcM ■ B 
indeed gives the average kinetic energy. The absolute minimum of the new curve, 



d(47iM ■ B) 
dH 



H 



0. 



(2) 



determines a critical field called H*. A consequence of the connection established here is 
that this field determines the maximum average kinetic energy, and from this one obtains 



the maximum root mean square current density, J ma x = \f(J 2 )- Thus using the present 
approach one obtains information about the maximum current density without the use of 
Bean's model. For this one needs the London expression F^ in « 27rA 2 (J 2 )/c 2 , where A is the 
penetration depth. Actually at H* the order parameter is not constant but we consider as 
such, to obtain the order of magnititude of this maximum current density circulating in the 
superconducting state. 

In this paper we also show that the ratio H*/H c2 becomes k independent for k larger than 
a few units. This property can be useful because it suggests that H c2 can be determined at 
a lower field, namely, H*. Frequently the upper critical field of many compounds cannot be 
reached experimentally just because it falls beyond the capabilities of a given experimental 
set up. 

The 4irM ■ B experimental curves were obtained here using a Niobium sample with an 
approximately spherical shape, mass m = 0.6487 g, and critical temperature T c = 8.5 K 
with a transition width of AT ~ 0.3K. The sample was obtained in an arc melt furnace from 
99.9 % pureness Nb wire. X-ray diffraction done in this sample shows the expected metallic 
Nb phase. Isothermal magnetization data was always obtained starting from a zero-field 
cooled procedure and using a commercial Quantum Design PPMS extraction magnetometer 
facility. 

New features and properties, not visible in the traditional way of plotting magnetization 
curves, are unveiled here by this 4irM -B plot. Some of them are explained in this paper, but 
there are others, such as the multi-minima structure in the irreversible region, that require 
further discussion and is not done here. In the reversible region, only one minimum remains 
and we find good agreement between theory and experiment. Close to T c the experimental 
and theoretical AirM ■ B curves show good agreement and the experimental ratio H*/H c2 is 
well predicted by the present theoretical models. Notice that the Virial Theorem, Eq. ([!]), 
holds for the microscopic theory ||, and for this reason the connection between the curve 
4nM ■ B and Fki n , remains valid beyond the GL theoretical framework of this paper. 

This paper is organized as follows. In Section || we show that the curve 4nM ■ B 
gives the average kinetic energy of the superconducting state and derive this curve both 
analytically and numerically in different theoretical situations. The field H* is calculated 
and, its properties, discussed here. In Section [HI] the experimental 4nM ■ B curves for the 
Nb sample are obtained, and the field H* determined. The phase diagram of this Nb sample 
is obtained with the curve H*(T) included. Finally we draw conclusions in Section [TV]. 
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II. THEORY 



Let us first obtain analytical expressions for the energy function AttM ■ B in two simple 
situations and determine there the critical field H*. The lower and upper critical fields are 
given by H c i = ($ /47rA 2 ) In A/£ and H c2 = $ /27r£ 2 , respectively, £ is the coherence length 
and k = A/£ the Ginzburg-Landau parameter. The first case investigated here is Abrikosov's 
theory which calculates the vortex lattice near H c2 in the framework of the GL free energy. 
There the magnetic induction is given by 

where (3 A = (|A| 4 )/(|A| 2 ) 2 is a parameter determined from the vortex lattice symmetry, A 
being the wavefunction that describes the superconducting state. One obtains that 

4vrM ■ B = - H f~ H (H - H f~ H a ), (4) 

is a parabola that vanishes at H c2 /(1 + (2k, 2 — l)(3 A ) and at H C 2. Although Abrikosov's 
theory is valid only near H c2l still, for k ^> y/2/2, the first vanishing field approaches 
2/(1 + 2/3^)$ /47tA 2 , which is of the same order of magnitude as the true H cl field. The 
absolute minimum of 4irM ■ B is achieved at the field 

H c2( 2f3 A K 2 -f3 A + 2 ^ 
H ~ — [ 2^-p A + l ] - (5) 

A simpler expression is obtained noting that f3 A is very close to one, which is true for 
the triangular lattice (f3 A = 1.1595953 ■• •): H* w (H c2 /2) * (1 + 1/k 2 ). Thus for any 
superconductor with k larger than a few units, this critical field is given by if ~ 0.5 H c2 . 

The second model is London theory, known to provide a good description of the super- 
conducting state in the intermediate field region far away from the critical fields. This is 
a fairly good approximation as long as variations of the order parameter near the vortex 
cores can be ignored, and so, the order parameter can be considered constant all over the 
superconductor. In this region the magnetic induction is very well approximated by |J 

B = H - H W«!EL. (6 ) 

inn 

Then the energy function becomes 

4nM.B = -BH cl ln{Hc2/B \ (7) 

In « 

whose minimum absolute is at the field value 

H* = ^-^±. (8) 
e In k 

Since H c i/H c2 = ln/t/2K 2 , notice, once more, that in case k is larger than a few units the 
last term can be neglected, rendering H* w 0.37H c2 . 
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The product 4irM ■ B is related to the average kinetic energy of the condensate, without 
any approximation, through the Virial Theorem (Eq. ([!])). 



4vrM • £? = -F km /2 - ((h - (h)Y) (9) 

is expressed here in reduced units, such that the thermodynamical fields and energies are 
scaled in the following way: h — > h/\/2H c , H — > H/y/2H c , and F kin — > F kin / (H 2 / 'Air) , H c 
being the thermodynamic field. The same scaling holds for the magnetic induction which is 
just the average local field, B = (h). These reduced units are convenient because to retrieve 
the above quantities at any given temperature T from a single temperature independent 
plot, it is enough to scale energies by H C (T) 2 /Arc, and fields by y/2 H C (T). 

To show that the function AnM-B yields the average kinetic energy, F kin , is just a matter 
of determining how the mean square deviation of the local field, ((h — (h)) 2 ), behaves as 
a function of k. Thus is done here for the GL theory using the iterative method of E. H. 
Brandt [|]]. This iterative procedure assumes an initial solution, from which one obtains 
very accurately, after a few iterative steps, the local magnetic field and the order parameter 
for any applied field in the mixed state. The AnM ■ B (continuous line) and F kin (points) 
curves, obtained using this method, are shown in Fig. [I] for several k, ranging from 0.75 to 
50. The conclusion is that these two curves, 4-7rM • B and F k i n /2 become very similar for 
large k. The field H* is obtained by determining the absolute minimum of the 4nM ■ B 
curves in Fig. [l|, and Fig. |2] summarizes these results as a H* /H c2 vs k plot. 

The difference between 4irM ■ B and F kin /2 is best quantified , as function of k, by 
two variance ratios introduced here and both computed at the field H*. We compute ((h — 
(h)) 2 ) /\4ttM ■ B\ H * for a few points («, ((h- (h)) 2 )/\ATrM ■ B\ H *): (0.75, 0.11), (0.85, 0.094), 
(1.0, 0.072), (1.5, 0.034), (2.0, 0.020), (3.0, 0.0088), (5.0, 0.0032), (10.0, 7.83 10~ 4 ), 
(20.0, 1.99 10- 4 ), (30.0, 8.8 10" 5 ), and (50.0, 3.2 10~ 5 ). Notice that for « > 3 the 
two curves can be regarded as the same with an error of less than one per cent in the 
neighborhood of H*. This error is not uniform for all fields, being smaller near H c2 and 
larger near H&. In Fig. ^|one has that H*/H c2 = 0.51 for k = 3.0, and essentially, for k > 3, 
H* is fairly well approximated by its limiting value of 0.5 H c2 . The inset of Fig. p shows the 
relative root mean square deviation of the local mag netic field, (((h - (h)) 2 )) 1/T / (h)\ h* , as 
a function of k. Again we see a sharp drop of this quantity as k increases. For instance for 
K = 3.0,(((h-(h)) 2 )) 1/2 /(h)\ H .= 0.0213. 

Analytical expressions for these two variances are obtained in the Abrikosov approxima- 
tion which starts from the well-known Abrikosov relation between the local field and the 
order parameter ||, h = H — 2tt\ A\ 2 /k, in reduced units. From it one obtains an interesting 
expression that relates thr variance to the magnetization: ((h — (h)) 2 ) = {AttM) 2 {(3a — 1)- 
From it one finds that at the field H*, given by Eq.(W), the following expressions: 



{{(h- {h)f)) l/2 (3 A -1 

771 \H* = { 7, 



,1/2 



(h) in v (3 A 4(3 a ^-4((3 a -1)k 2 + ((3 a ~1)' 

and, 

((h-(h)) 2 )t (3 A -1 



(10) 



\4ttM-B\ W ' 2(3 a k 2 + 1-(3 a 



(11) 
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Both expressions display a 1/k 2 behavior for large k. For the triangular lattice, one 
obtains that (((h - (h)) 2 )) 1/2 / (h) \ H . w 0.172256/k 2 , and ((h - (h)) 2 )/\4irM ■ B\ H , w 
0.068815/k 2 , which agree quite well with the numerical method results. We conclude that 
the curve 4nM ■ B can be obtained from the average kinetic energy for k > k c , with a 
precision better than one percent for k c pa 3. Next we study a Niobium sample and obtain 
the curve 4ttM ■ B experimentally, starting from the isothermal magnetization curves. 



III. EXPERIMENT 

To exemplify the present ideas we have studied superconductor Niobium, which along 
the years has been an important reference for the understanding of type II superconducting 
properties . The observation of the irreversibility line in pure Nb |7j] as observed in high-T c 
superconductors H], has renewed the interest on this superconductor. Also, effects which are 
easier to study and understand in Nb, such as surface pinning ||10|| , vortex avalanches [11 



and peak effect |12|| , has shed some light on the understanding of similar effects in high-Tc 
superconductors. 

The sample has critical temperature below of pure Nb (T c = 9.2K), indicating that 
it contains non-magnetic impurities. This causes no problem here, on the contrary, the 
impurities turn the compound into a more type II superconductor, with a larger value of k. 
Very pure Nb has K pa 1.0, whereas we find here k pa 4 for this sample, when comparing 
the present theory to the experiment. 

Next we discuss our data as shown in the figures. 

The transition properties of this sample are seen in Fig. |3|, which shows a rather small 
reversible temperature region. For fields below 4.0 kOe, the magnetization becomes essen- 
tially irreversible. Values of H c2 (T) are extracted from this diagram following the standard 
linear extrapolation of the magnetization curve. 

Fig. ^ shows a total of nine isothermal magnetization cycles, 0.5 K apart from each 
other, for temperatures ranging from 3.5 K to 8 K. Notice that in the neighborhood of zero 
applied field all the magnetization curves overlap forming one single temperature indepen- 
dent line. The slope of this line defines the demagnetization factor. For finite size samples 
the demagnetization factor must be included []14|| ; for a sphere this factor is 1/3 such that 
B = H in + 47rM/3, where H in — H + 4irM/3 is the internal field. Our sample is an ellipsoid 
with axis 2n = 4.7 mm and 2r 2 = 5.3 mm, whose sphere with equivalent volume has radius 
5.1 mm. We have experimentally determined the demagnetization of our sample using the 
following procedure. We assume that B = H + 2dM where d is a factor to be determined 
from our data such that, in the Meissner phase, the magnetic induction must vanish. We 
find that, in this region, a linear extrapolation of the magnetization data for all AwM vs H 
curves produces the same value of 2d = 4.16, independent of temperature. 

Another important issue here is the removal of the background magnetization as a func- 
tion of H for all temperatures. The construction of the curve AirM ■ B is only useful, e.g. to 
determine the critical field H*, if all magnetic contributions that are not of superconducting 
origin are totally removed from the original isothermal M vs H data. For the Nb sample 
the following procedure was applied. The raw magnetization data displays a linear non-zero 
positive slope as a function of H that goes beyond the H c2 field. This curve is universal and 
temperature independent. We subtract from the raw magnetization data at all temperatures 
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this single linear curve M vs H. This is possible because the background magnetization is 
temperature independent, and just a single measurement of the paramagnetic background 
signal close or above T c is enough to subtract the unwanted background for any tempera- 
ture. It is well-known that this source of magnetic signal is the Pauli paramagnetism of the 
normal electrons. 

Selected isothermal curves 2dMB vs H are shown for three different temperatures, cor- 
responding to T = 4.5, 6, and 8 K (Figs. |5], ^ and |7|, respectively). The basic features of 
the M vs H curves, such as the critical fields and the peak effect, are also present in the 
2dMB vs H curves, although in a somewhat different way. However the 2dMB shows a 
new feature, not noticed in the traditional M vs H plot, namely the new minimum at H* . 
For T = 4.5 K the absolute minimum of the curve 2dMB occurs at H = 6.0 kOe, a field 
associated with the peak effect. Still the H* is observable there. However for increasing 
temperature, T = 6 K, H* becomes the absolute minimum and the peak effect turns into a 
local minimum, in this case around 4 kOe. Finally at T = 8K only one minimum remains, 
H*, and the peak effect is essentially absent. Notice that two H* fields are defined, measured 
in ascending (zero field cooled - zfc) and descending (field cooled - fc) fields. The zfc and fc 
H* fields are separated by approximately 500 Oe for the T = 4.5 K curve whereas for 8 K 
they become the same. 

In Fig. [8] we compare the k — 3, 4, and 5 theoretical and the T = 8.0 K experimental 
curves obtained in this paper. All 2dMB data sets are normalized along the 2dMB axis to 
their maximum absolute values, achieved at H*. Along the H axis, the theoretical curves 
are scaled by the inverse of their upper critical fields, such as in Fig. [l|. However for the 
experimental curve we adopted another normalization procedure, namely we scaled fields 
by 796.98 H due to the following reason. Our goal is to shift the 2dMB experimental 
curve along the H axis in search of the best theoretical fit, which we find to occur when 
experimental and theoretical H* fields coincide. The experimental absolute maximum of 
\2dMB\ is reached at H* = 329.81 Oe. Next we apply the scale transformation aH/329.81 
and search for a that gives the best agreement with the theoretical curves. This happens at 
a = 0.46, very close to the theoretical expected ratio of a = H*/H c2 ~ 0.5 of the theoretical 
curve, described in Fig. ^|. Notice that the three theoretical curves coincide near H c2 , as 
they should, because of the normalization, but not near H c \. From this figure we conclude 
that the theoretical k = 4 curve provides the best fitting to the experimental data. It is 
worth mentioning that this value of k agrees with its experimental estimate derived from 
the ratio of the critical fields. H c2 /H c i. 

Fig. |9] presents the phase diagram, obtained from the values of H c2 that follows from the 
curves of Fig. Also plotted in this figure are the values of H* found in the increasing and 
decreasing field branches of plots, such as shown in Figs. ||, ||, and |7|. 

We find a quite perfect match between the experimental values of H* and H c2 (T) with 
the expression obtained from the theories previously discussed. Thus for this Nb sample we 
find that the curve H c2 /2 fits fairly well the experimental H* curve. 

We estimate the maximum root mean square average current, J m , for T = 8 K. Firstly 
one must obtain the London penetration depth from the lower critical field for this tempera- 
ture: H c \ ~ 18 Oe. We obtain that A ~ 0.36 fim and taking that the maximum of \2dM ■ B\ 
is reached at 3.88 10 3 G 2 , it follows that J m ~ 4.2 10 5 A. cm -2 . Notice that this maximum 
mean square average current density is being obtained in the reversible regime. 
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In summary we find here that this AirM ■ B curve naturally enhances many aspects of 
traditional 4ttM vs H plots, such as the peak effect, the first entry field (H c i), which is also 
much more clearly resolved in this kind of plot. Obviously the theories considered in the 
previous section do not describe the many features of the experimental curves 2dMB vs H, 
especially when irreversible processes are present. This is expected because those theories 
do not take into account pinning effects. When pinning effects become small and there is 
reversibility in the magnetization, we found agreement between theory and experiment. It 
is always possible to cast isothermal magnetic data as a AtcM ■ B curve, but the question 
is whether the present connection between the 4nM ■ B curve and F^ in still holds when 
irreversible effects are present, such as the peak effect. We believe that this connection indeed 
remains valid in the irreversible region in some large k limit. Obviously this connection 
should always be regarded approximately and corrections are expected due to pinning. We 
believe that these corrections should be small, based on the fact that the Virial Theorem 
gets little, or no modification at all || when some types of defects H JT3|JT^[] are present. For 
temperatures close to T c , where irreversible effects become small, we found that the present 
ideas hold and the average kinetic energy interpretation of the curve AuM ■ B is applicable. 



IV. CONCLUSIONS 

In conclusion the present proposal of casting isothermal magnetization data as the prod- 
uct of the magnetization times the magnetic induction is useful because it directly gives the 
average kinetic energy of the superconducting state for any applied field. We have stud- 
ied here this curve both theoretically and experimentally. The field of maximum average 
kinetic energy in the condensate is just the minimum point of the new curve. From this 
field we estimate the maximum mean square average current density. We show using the 
Ginzburg Landau theory that this field value is half the upper critical field value for large 
k superconductors. Thus its experimental determination can help to understand proper- 
ties of the upper critical field which often falls beyond the experimental capabilities of a 
given equipment. We found analytical expressions for this curve in two special limits of 
the Ginzburg-Landau theory, namely the Abrikosov and the London limits. In the general 
situation we have obtained this curve numerically for several values of the Ginzburg-Landau 
constant, using an efficient iterative method [|J to solve the Ginzburg-Landau theory. We 
find that for a Ginzburg-Landau constant larger than a few units the mean square deviation 
of the local magnetic field becomes extremely small. Only for such parameter values the 
product magnetization times the magnetic induction gives the average kinetic energy. 

In summary we find the present proposal of plotting the isothermal magnetization useful 
because it provides a new insight into the properties of the mixed state. 
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FIGURES 



FIG. 1. Two theoretical curves, AirM ■ B (continuous) and Ffcj„/2 (dots), versus H are shown 
here for several k, namely, 0.75, 0.85, 1.0, 1.2, 1.5, 2.0, 3.0, and 5.0. For higher k, 10, 20, 30, and 
50, the two curves, AttM ■ B (continuous) and (points) are shown in the inset. The field H*, 

the absolute minimum of the 4irM • B curve, is suggested for a few curves, as examples. The inset 
clearly shows that for large n H*/H C 2 ~ 0.5. 

FIG. 2. The field H*, the absolute minimum of the 4itM ■ B curve, is plotted here as a function 
of k. For k larger than a few units the ratio H*/H C 2 approaches 0.5. The inset shows that the 
root mean square deviation of the local field vanishes for k larger than a few units. 



FIG. 3. The zfc and fc magnetization versus temperature curves are shown here for several 
applied fields, ranging from 1.0 kOe to 9.0 kOe every 1.0 kOe. Below 4.0 kOe the curves can be 
considered as reversibles. 



FIG. 4. Isothermal magnetization cycles measured in ascending (ZFC) and descending (FC) 
fields are shown, starting from 3.5 K to 8 K every 0.5 K. The overlap of the inital part of the cycles 
into a single straight line defines the demagnetization factor of the sample. 



FIG. 5. A selected isothermal magnetization cycle is plotted here both in the traditional and 
in the new ways, for the temperature of 4.5 K. The field H* is hardly noticeable whereas the peak 
effect is the absolute minimum of the 2dMB vs H curve. 



FIG. 6. A selected isothermal magnetization cycle is plotted here both in the traditional and 
in the new ways, for the temperature of 6.0 K. The field H* is the absolute minimum, whereas the 
peak effect is just a local one, for the IdMB vs H curve. 



FIG. 7. A selected isothermal magnetization cycle is plotted here both in the traditional and 
in the new ways, for the temperature of 8.0 K. The field H* is the absolute minimum of the IdMB 
vs H curve, which is essentially totally reversible. 
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FIG. 8. This plot shows the theoretical curves (n = 3, 4 and 5) and that, among them, the 
k = 4 curve provides the best fitting to the T = 8K experimental curve. All four curves were 
normalized to their maximum absoulte value. Along the H axis the three theoretical curves were 
scaled by their H C 2, and the experimental curve was scaled such that its H* field is 0.46 away from 
the theoretical normalized H C 2 field. 



FIG. 9. The phase diagram of the Nb sample, containing its several critical fields, is shown here. 
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